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consideration  in  the  selection  of  specific  systems  for  study,  particularly 
multilayer  integrated  circuit  applications. 

UO  A 

In  this  period,  an  in-situ  argon  ion-gun  cleaning  technique  has  been  used 
to  strip  the  native  SiO?  film  from  the  Si  substrates  prior  to  the  actual  depo- 
sition  of  Si^and  SiC  films.  Etch  rate  data  for  the  native  oxide  on  SiO^  and 
the  relative  film  thickness  for  reactively  sputtered  SiO-  films  as  a  function 

f*)i  * 

of  substrate  temperature  for  etched  substrates, ha*  been  gathered.  The  relative 
deposition  rate  and  film  texture  for  single-source  sputtered  SiC  films  as  a 
function  of  substrate  type,  orientation  and  film  texture  was  obtained.  Explor¬ 
atory  films  of  SiQ  .  W0„  and  NbCn  QQ  have  been  fabricated  using  the  sequential 
deposition  of  atomic  layers  technique. 

In  the  theoretical  area,  the  CNDO/2  technique  has  been  used  to  determine 
the  energies  of  formation  for  bridging  oxygen  vacancies*.  V«,  and  Frenkel  de- 

•r  ^ 

fects,  Vg  +  Oj,  for  a  pair  of  silica  structural  units  (dimer  and  tetramer 
ring).  In  addition,  calculations  have  been  made  for  the  adsorption  of  dif¬ 
ferent  molecular  chemical  species  (in  the  parallel  configuration)  at  the 
bridging  oxygen  site  of  a  silica  dimer.*  In  the  crystallographic  modeling  of 
silica  structures,  energetic  calculations  have  indicated  that  the  0-0  distance 
provides  an  important  restrictive  criterion  for  the  configurational  changes  of 
the  Si 0^  tetrahedra  in  vitreous  silica.  Using  this  criterion  as  a  selective 
sieve,  allowed  bridge  bond  angles  for  the  trimer  unit  were  determined  for  a 
pool  of  1200  randomly  selected  examples.  Only  61  examples  passed  the  sieve 
test  and  these  had  a  distribution  of  bridge  bond  angles  centered  at  ~  145° 
with  a  half  width  of -30°. 
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INTRODUCTION 


During  this  fourth  six-month  period  of  the  program,  attention  was 
given  to  stripping  the  native  oxide  from  the  Si  substrates  prior  to  Si 0X 
and  SiC  deposition.  The  SiC  deposition  was  from  a  single  source  and  a 
beginning  excursion  was  made  into  the  use  of  the  layer  by  layer  method  of 
film  formation  for  two  oxides.  The  theoretical  emphasis  continued  to 
focus  on  both  the  structural  and  the  energetic  aspects  of  SiOx. 
Consideration  of  the  trimer  unit  appears  to  be  the  key  needed  for 
understanding  vitreous  Si O2  formation  on  a  macroscale.  On  a  microscale, 
insights  into  bridging  oxygen  vacancy,  Vq,  and  Frenkel  defect,  Vq+Oj, 
formation  energetics  have  been  gathered.  To  investigate  the  electrical 
properties  of  these  films,  attention  has  been  given  to  the  development  of 
a  new  data  acquisition  system  based  upon  the  ac  electrical  conductance 
method. 

A.  Experimental  Results 
1 .  Synthesis  of  SiOv 

Thus  far  in  this  program,  synthesis  experiments  via  reactive 
sputtering  have  consisted  of  continuous  process  experiments,  i.e.,  both  O2 
and  Si  were  incident  on  nominally  room  temperature  substrates  concurrently 
with  reaction  taking  place  to  form  SiOx  films.  The  basic  result  of  this 
work  is  that  films  of  uniform  composition  can  be  deposited  with  control  of 
stoichiometry  in  the  range  of  0.3  <  x  <  2  .  During  a  given  synthesis 
run,  the  substrate  temperature  did  not  float  higher  than  50°C. 

During  the  recent  period,  work  was  performed  to  extend  and  refine 
these  experiments  so  as  to  increase  our  understanding  of  the  process  and 


of  the  types  of  materials  which  can  be  synthesized.  Two  different 
approaches  were  taken  to  the  film  deposition,  the  first  following  directly 
on  the  earlier  work.  In  the  first  approach,  the  system  was  improved  by 
adding  an  in-situ  cleaning  technique  using  an  argon  ion-gun  for  stripping 
the  native  Si O2  film  from  the  substrate  prior  to  the  actual  deposition  of 
the  reactively  sputtered  Si 0X  films.  In  addition, a  new  O2  gas  dispersion 
ring  of  rectangular  shape  was  developed  for  processing  rectangular 
substrates.  These  new  components  are  illustrated  in  the  chamber  schematic 
of  Fig.  1.  Some  ion-etching  data  is  given  in  Table  I  where  we  see  that  a 
substantial  etching  -ate  occurs  for  a  2.0  kV  anode  voltage  located  4"  from 
the  substrate  with  a  4  yin  argon  pressure  and  that  the  procedure  is  very 
sensitive  to  both  voltage  and  distance.  Deposition  experiments  on  Si (111) 
substrates,  ion  etched  at  750  V  and  2  mA  for  40  minutes,  were  carried  out 
on  heated  substrates  with  the  film  thickness  results  given  in  Table  II. 
There  is  ',on;e  question  concerning  the  complete  removal  of  the  native  oxide 
for  these  experiments.  Using  alpha-step  and  ellipsometry  measurements, 
the  relative  film  thickness  was  measured  for  these  films  and  we  see  from 
Table  II  that  the  film  thickness  decreases  by  a  factor  of  2  between 
deposition  at  25°C  and  660°C.  We  note  also  that  the  film  stoichiometry 
becomes  increasingly  Si-rich  as  the  substrate  temperature  is  increased. 
This  may  be  the  reason  why  the  plot  of  film  thickness  vs.  1/T  in  Fig.  2 
does  not  reveal  a  linear  relationship. 

An  additional  series  of  experiments  was  performed  wherein  the  Si  and 
O2  were  not  concurrently  incident  on  the  substrate.  The  chamber  system 
was  modified  so  that  the  substrate  was  sequentially  exposed  to  Si  and  to 
O2  so  that,  in  principle,  a  monolayer  of  amorphous  Si  was  first  deposited 
and  then  subsequently  oxidized  by  exposure  to  an  O2  atmosphere  for  a  time 
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needed  to  produce  a  fraction  of  a  monolayer  of  O2  and  then  the  entire 
process  repeated  many  times  to  produce  an  Si 0X  film.  Some  experimental 
data  with  this  method  is  presented  in  Table  III.  More  will  be  said  about 
this  general  technique  in  the  next  section. 

2.  Synthesis  of  Compounds  by  Atomic  Layer  Growth 

Over  the  past  decade,  a  new  technology  for  the  growth  of  high 
perfection  crystalline  materials  has  been  developed.  This  technique, 
"Molecular  Beam  Epitaxy,"  allows  the  synthesis  of  a  wide  variety  of  unique 
structural  deposits.  It  has  characteristically  been  applied  to  GaAs  and 
to  (Ga.Az)As  .  In  this  process,  an  overpressure  of  As  is  maintained  at 
the  substrate  and  the  growth  of  the  deposit  is  defined  by  the  arrival  rate 
of  Ga  and  At  .  Atomic  layer  growth  is  observed  and  smooth  interfaces  can 
be  formed  which  are  abrupt  on  an  atomic  or  near  atomic  scale. 

In  our  laboratory,  a  technique  using  sputtering  technology  has  been 
developed  that  allows  sequential  deposition  of  atomic  species  that  can  be 
controlled  on  a  monolayer  level  and  which  is  uniform  for  many  hundreds  of 
layers.  In  this  process,  either  (1)  separate  sputter  deposition  sources 
are  used  for  each  constituent  or  (2)  a  separate  region  is  provided  for  the 
substrate  to  be  exposed  to  reactive  gases.  The  samples  are  mounted  on  a 
turntable  rotating  at  a  fixed  rate  controlled  to  0.1%  and  pass  beneath 
sputter  sources  so  that  sequentially  deposited  layers  of  the  various 
constituents  are  incident  on  the  substrate.  Individual  layer  thickness  is 
controlled  by  the  rotation  of  the  substrate  table  and  the  deposition  rates 
of  the  various  sources  used.  In  the  case  of  a  reactive  gas,  it  appears 
that  the  partial  pressure  of  the  gas  over  the  substrate  is  the  most 
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important  factor  and,  when  coupled  with  the  rate  of  substrate  motion, 
determines  the  exposure  of  the  deposition  surface. 

At  this  time,  exploratory  experiments  have  been  performed  on  the 
synthesis  of  SiOx  and  W0X  films.  It  has  been  clearly  demonstrated  that 
large  area  films  of  uniform  thickness  and  stoichiometry  can  be  deposited 
by  this  technique.  In  addition,  it  has  been  demonstrated  that  NbCggg  (of 
the  NaC*  structure)  can  be  synthesized  onto  room  temperature  substrates 
and  exhibits  the  superconducting  transition  temperature  expected  for  bulk 
material  of  this  composition. 

3.  Synthesis  of  SiCvThin  Films 

In  the  last  report^)  we  described  our  work  on  two-source 
deposition  of  SiC  from  Si  and  C  targets.  This  type  of  deposition 
necessarily  results  in  a  film  of  varying  concentration  from  C-rich  at  one 
end  of  the  substrate  to  Si -rich  at  the  other.  However,  we  are  primarily 
interested  in  stoichiometric  3-SiC  and  ways  of  improving  the 
crystallinity  of  sputtered  SiC  films.  Thus,  we  have  changed  from  the  two- 
source  to  a  single  source  arrangement.  By  using  a  disk  of  pressed 
o-SiC  particles  as  the  target,  we  are  able  to  sputter  deposit 
stoichiometric  B-SiC  thin  films. 

The  actual  deposition  process  has  been  improved  by  adding  an  in-situ 
cleaning  technique.  The  new  deposition  configuration  shown  in  Fig.  3 
includes  the  SiC  source,  the  argon-ion  gun  for  surface  cleaning  and  a  tin 
source  which  can  be  moved  to  any  desired  location  in  the  chamber.  The  top 
view  of  the  system  shows  the  substrate  travel  path  for  a  deposition 
cycle.  After  pumpdown  to  about  4  x  10"®  Torr,  the  substrate  heater  is 
slowly  ramped  to  the  deposition  temperature.  The  substrate  is  rotated 
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under  the  argon-ion  gun  and  surface  sputter  cleaning  is  initiated.  It  is 
essential  to  prepare  a  clean  surface  for  good  SiC  epitaxy  to  occur  so  the 
native  SiC^  on  the  Si  substrate  must  be  completely  removed.  This  involves 
up  to  45  minutes  of  argon-ion  bombardment  at  2  mA  and  750  V,  a  voltage 
considered  to  be  low  enough  to  prevent  significant  surface  damage  by  ion 
bombardment.  Before  being  loaded  into  the  chamber,  the  Si  substrates  were 
put  through  the  industry  standard  "RCA"  clean  treatment.  Near  the  end  of 
the  sputtering  cycle,  the  SiC  source  is  activated  and  the  substrate 
quickly  moved  under  the  activated  SiC  target.  The  tin  source  will  be 
placed  so  that  the  Sn  is  co-sputtered  onto  the  substrate  surface  along 
with  the  SiC,  but  it  will  be  arranged  so  that  the  substrate  is  rotated 
through  the  tin  just  before  being  located  under  the  SiC  target. 


Deposition  Rate 


Not  only  is  a  knowledge  of  the  deposition  rate  important  in  order  to 
generate  samples  of  a  desired  thickness  but  it  also  can  be  used  to  compare 
alternate  deposition  techniques.  The  deposition  rate  can  even  give 
insight  into  the  degree  of  crystallinity  of  the  deposited  film.  For  the 
dual  source  method  at  440  W  and  a  substrate  temperature  of  700°C,  the 
deposition  rate  was  ~  2  A/sec  .  The  single-source  deposition  is  almost 
twice  as  fast  (-  4  A/sec)  .  The  following  trends  have  been  observed  for 
the  deposition  rate,  D.R. 


D.R.  (Rm  T,  A^C^)  ■  1.2  D.R.  (elev  T,  Aj^Oj)  > 

D.R.  (Rm  T,  Si (111  or  100))  ~  D.R.  (elev  T,  Si (111  or  100)), 

D.R.  (Sim)  a  1.2  D.R.  (Sijgg), 


D.R.  (Si)  »  1.4  D.R.  (A22O3)  at  Rm  T  , 

D.R.  (Si)  «  2.0  0.R .  (Ai203)  at  700°C  . 

Any  temperature  or  substrate  effects  might  be  attributed  to  differences  in 
film  crystallinity.  Amorphous  films  are  expected  at  low  temperatures  and 
mismatched  substrates,  whereas  crystalline  fi lms  are  deposited  above  600°C 
and  on  well  matched  substrates.  The  orientation  effect  indicates  the 
higher  compatibility  with  the  (111)  substrate  and  this  appears 
dramatically  in  the  x-ray  analysis.  The  substrate  type  effect  illustrates 
the  much  greater  affinity  of  the  SiC  for  Si  than  for  even  though 

there  is  a  smaller  disregistry  on  the  (lT02)  A^O-j  than  on  the  Si 
substrates.  The  greatest  deposition  rate  is  observed  on  the  substrate 
surface  having  the  same  symmetry  as  the  ''natural"  or  "preferred" 
orientation  of  the  deposited  B-SiC  . 

X-Ray  Diffractometry 

Our  primary  analytical  technique  to  determine  the  perfection  of  the 
sputtered  films  is  X-ray  diffractometry.  Me  use  a  Read  camera  to  evaluate 
the  extent  of  preferred  orientation.  As  with  the  dual  source  samples,  a 
strongly  preferred  orientation  of  the  single  source  films  was  found.  Only 
the  (111)  B-SiC  reflection  is  seen.  A  comparison  of  one  deposition 
(#178),  which  had  a  specular  surface  and  gave  the  strongest  X-ray  peak, 
with  other  depositions  is  instructive.  The  data  are  presented  in  Table 
IV.  The  relative  peak  heights  were  derived  from  X-ray  diffraction  peak 
heights  accounting  for  peak  broadening  (same  in  all  cases )  and  SiC  film 
thickness.  Ignoring  the  difference  in  the  3“  and  4"  variable  which  is 
considered  to  be  a  second  order  effect,  the  data  show  that  deposition  on 
clean  (111)  Si  greatly  enhances  film  crystallinity  for  the  source  operated 
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at  600  W.  Read  camera  photographs  showed  a  greatly  enhanced  preferred 
orientation  for  sample  178.  In  addition,  the  sample's  edges  and  surface 
were  featureless  under  optical  and  scanning  electron  microscopy.  By 
comparison,  the  more  typical  film  surface  appears  rough  and  scanning 
electron  microscopy  of  the  dual  source  samples  revealed  columnar  growth 
with  a  grain  diameter  -  100  A  . 

Electron  Microprobe 

Figure  4  presents  electron  microprobe  results  for  three  single  source 
samples  and  one  dual  source  sample.  These  samples  were  analyzed 
sequentially  and  were  referenced  to  the  same  SiC  microprobe  standard, 
which  is  a  polished  piece  of  a-SiC  .  The  films  are  0-Si C  deposited  on 
SOS  quality  sapphire  at  Rm  T,  700°C  and  820°C.  The  single  source 
deposition  at  Rm  T  and  700°C  appear  to  be  stoichiometric  SiC  and  uniform, 
while  the  820°C  sample  seems  to  be  off-stoichiometry  on  the  C-rich  side 
but  fairly  uniform  in  composition  across  the  sample. 

Refractive  Index  and  Optical  Transmission 

Refractive  index  (R.I.)  measu, ements  have  been  made  on  a  Gaertner 
L 1 1 7  ellipsometer  equipped  with  a  HeNe  laser.  Films  of  SiC  deposited  at 
600  W  appear  to  have  R.I.  ~  3.2  while  films  deposited  at  200  W  have 
R.I.  ~  2.6  (literature  values  give  R.I.  ~  2.654  -  2.697  for  s-SiC)  . 

The  R.I.  measurements  on  the  600  W  samples  are  somewhat  suspect  because 
the  ellipsometer  null  was  quite  broad. 

We  have  been  hoping  to  obtain  definitive  measurements  of  the  band  gap 
by  optical  transmission.  The  sample  evaluated  previously  as  having 
~  3.0  eV  band  gap  was  re-evaluated  on  another  instrument  (Beckman  DK-2A 
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Spectrophotometer)  and  predicted  to  have  a  value  of  2.5  when  one  takes 
account  of  film  thickness,  transmission  to  absorbance  conversion  and 
assuming  a  direct  band  gap  (assumption  of  an  indirect  gap  gives  a  very 
low-value).  The  literature  values  vary  with  recent  reports  citing 
-  3.4  for  3-Si  and  -  3  for  a-SiC  . 

In  this  measurement  technique,  we  are  concerned  with  substrate 
scattering  since,  to  date,  our  optical  transmission  substrates  have  been 
only  singly  polished  sapphire.  Doubly  polished  sapphire  has  been  on  order 
and  is  expected  soon.  We  hope  to  obtain  better  results  with  these  new 
substrates. 

B.  Theoretical  Results 

1 .  ac  Electrical  Conductance  Measurement  Techniques 

In  the  fabrication  of  single  layer  or  multilayer  structures 
for  device  applications,  ultimate  characterization  is  usually  that  of 
electrical  properties.  However,  two  primary  inhomogeneities  that  are  not 
effectively  discriminated  by  the  present  dc  measurement  techniques  are 
those  at  the  interfaces  in  the  system  and  those  between  the  electrodes  and 
the  sample.  Here,  ac  measurements  of  the  dielectric  response  of  the 
film/substrate  system  can  be  invaluable  for  detecting,  and  separating  out, 
various  contributions  to  various  pol ari zati ons  in  the  system  under 
study.  In  spite  of  the  vast  amount  of  studies  on  Si-0  and  other  film 
systems,  the  significance  of  ac  conductivity  techniques  has  not  been 
widely  appreciated,  primarily  because  of  the  difficulty  in  translating 
frequency -dependent  behavior  into  actual  physical  mechanisms. 

This  is  because  the  past  techniques  have  been  either  (1)  accurate  but 
slow  or  (2)  fast  but  inaccurate  and  because  no  generalized  equivalent 
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circuit  model  has  existed  via  which  various  contributions  to  the 


dielectric  response  spectrum  could  be  assigned  to  specific  circuit 
parameters. 

From  our  studies  in  this  area,  we  have  devised  a  system  which 
promises  to  overcome  all  of  the  above  limitations.  It  allows  rapid 
collection,  storage  and  analysis  of  large  amounts  of  data  at  the  site  of 
the  film  experiment.  The  device  is  called  a  dielectric  response 
analyzer.  It  will  measure  comparative  gain  and 

pnase  angle  between  input  and  output  sinusoidal  stimuli  to  the  system  by 
means  of  a  dual  channel  recorder.  The  rapid  sampling  features  of  the 
recorder  will  extend  the  upper  frequency  range  to- 5  MHz  and  the  lower 
limit  to  -  10‘3  Hz. 

In  addition  to  providing  a  rapid  and  accurate  measuring  system  for 
the  following  areas:  (1)  Research  on  (a)  annealing  rates,  (b)  photo¬ 
stimulation  response,  (c)  voltage  bias  effects,  (d)  charge  carrier 
species,  (e)  diffusion  coefficients  and  (f)  electrode  effects,  plus 
(2)  Quality  Control  via  (a)  conductance,  (b)  capacitance,  (c)  dielectric 
constant,  (d)  interface  integrity  and  (e)  microstructure,  the  new 
technique  will  allow  us  to  test  the  feasibility  of  non-contacting 
electrodes  (since  the  dielectric  response  is  usually  a  measure  of 
resonance  rather  than  of  macroscopic  current  flow).  Such  a  system  would 
allow  the  possibility  of  non-invasively  monitoring  electrical  changes  in 
films  under  a  very  broad  range  of  environmental  conditions. 

2.  Modeling  of  Silica  Structures 

In  our  last  report  for  this  topic^),  our  goal  was  to 
identify  key  parameters  for  describing  the  spatial  structure  of  the  known 
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polymorphs  of  silica  without  explicitly  invoking  the  unit  cell  dimensions 
and  space  groups.  Instead,  it  was  in  terms  of  Si O2  tetrahedra,  the 
angular  parameters  describing  the  orientation  of  adjacent  tetrahedra  and 
of  closed  rings  made  up  of  bridging  oxygens.  Procedures  have  been 
developed  for  building  ring  structures  of  different  size  and  shape;  i.e., 
finding  the  geometrical  constraints  in  the  dimer  parameters  that  would 
lead  to  closed  rings.  We  have  sought  additional  restrictive  criteria  on 
the  parameters  involved  in  building  ring  structures  and  have  focused  on 
the  intermolecular  potential  energy  changes  associated  with  changing  0-0 
and  Si-0  distances.  Trends  in  the  energy  variation  of  the  resulting 
structure  as  a  function  of  the  parameters  have  been  studied.  Attention 
has  been  focused  mainly  on  the  coulombic  energy  contribution  to  date; 
however,  semi -empirical  procedures  for  including  covalent  energy  changes 
and  repulsive  energy  changes  are  also  under  development. 

The  major  difficulty  in  this  type  of  quantitative  assessment  is  in 
finding  a  theoretical  approach  that  is  simple  enough  to  be  manageable  and 
yet  accurate  enough  to  give  meaningful  insight  to  the  problem.  In  the 
present  case,  the  dilemma  of  the  simple  approach  is  “what  is  the 
partitioning  of  the  total  energy,  Ej,  into  the  ionic  part,  E-j ,  and  the 
covalent  part,  Ec,  and  what  are  the  proper  partial  charges  to  be  used 
in  calculating  the  ionic  part?".  One  such  simple  approach  is  given  by 
Sanderson^)  who  was  very  successful  in  using  it  to  estimate  the 
atomization  energy  of  the  polymer  (Si O2 ) n  •  He  defined  blending 
coefficients  T^  (ionic)  and  Tc  (covalent)  such  that  T,-  +  Tc  =  1  with  a 
total  energy  of 


He  also  defined  partial  charges  D( Si )  =  +0.40  and  D(0)  =  -0.20,  for  Si  and 
0,  respectively,  and  defined  Tj  as 


T.  =  \  CD(Si)  -  0(0)] 


(2) 


Assuming  an  Si-0  distance  of  1.61  A,  this  led  to  =  0.3  and  Tc  *  0.7 
with  an  Si-0  single  order  ionic  bond  energy  of  E^  ■  61.9  kcal  per  mole 
(KPM)  and  a  Si-0  single  order  covalent  bond  energy  of  49.3  KPM.  This 

yields  an  Ey  of  111.2  KPM  per  SiO  bond  or  an  atomization  energy  for  the 

Si O2  polymer  of  444.8  KPM,  which  is  in  excellent  agreement  with  the 

experimental  value  of  445.8  KPM.  This  all  seems  very  fine  except  there 

is  an  appreciable  body  of  data  pointing  to  the  charge  on  the  oxygen  in 
Si 02  being  ~  -0.5  to  -0.8^3,4^  . 

Focusing  our  attention  solely  on  the  coulomb  energy,  CE,  it  is 
expressed  in  KPM  for  an  ion  pair  by 


CE 


332  Q1Q2 
R12K 


(3) 


where  332  is  the  conversion  factor  to  KPM  if  the  charges  on  the  ions 
and  Q2  are  in  electron  units,  the  separation  between  the  ions,  Rj2  is  in 
angstrom  units  and  K  is  the  effective  dielectric  constant.  The  parameter 
K  just  factors  in  the  polarization  effect  of  the  macroscopic  environment 
surrounding  ion  #2  due  to  the  electric  field  of  ion  #1.  We  shall  follow 
this  procedure  to  estimate  CE  for  a  variety  of  silica  clusters;  (i)  an 
isolated  tetrahedron  with  no  bridging  oxygens,  (ii)  an  end  unit 
tetrahedron  with  one  bridging  oxygen;  i.e.,  an  isolated  dimer,  (iii)  an 
internal  tetrahedral  unit  in  a  closed  ring  with  two  bridging  oxygens; 


i.e.,  an  isolated  ring  and  (iv)  a  tetrahedral  unit  in  an  extended  volume 
with  four  bridging  oxygens.  Since  we  are  comparing  isolated  species  with 
no  macroscopic  silica  environment,  it  is  reasonable  to  set  K  =*  1  .  This 
will  give  us  the  relative  coulomb  energy  for  a  tetrahedron  in  all 
possible  structural  configurations  for  the  silica  network.  We  can  ensure 
charge  neutrality  and  stoichiometry  of  the  Si O2  provided  we  assume  a 
weighting  of  1  for  all  Si,  a  weighting  of  1  for  all  bridging  0  and  a 
weighting  of  1/2  for  all  non-bridging  0. 

In  the  calculation  to  follow,  we  assume  a  regular  tetrahedron  (TET) 
of  0  with  a  Si  interstitial  and  an  Si-0  distance  of  1.60  A  which  yields 
an  0-0  distance  of  2.16  A  .  We  assume  a  negative  charge  on  the  oxygen 
in  electron  units  of  Q2  and  a  positive  charge  on  the  Si  of  +  2Q2.  We 
find  the  total  CE  of  the  unit  to  be  given  by 

CE  =  332Q2  )  2  F  S+  4£  B  +  £  ft]  KPM  (4a) 

(  SiO  SiSi  00  ) 

where 


0  =  WTX  •  WT2/R12  (4b) 

Considering  only  a  single  tetrahedron  (no  interactions  with  neighbors), 
we  have  CE/Q2  =  -639,  -752  and  -832  KPM  for  an  isolated  unit  (no  BRO), 

an  end  unit  (1  BRO)  and  an  internal  unit  (2  BRO),  respectively.  The  end 
unit  would  correspond  to  a  dimer  unit  while  the  internal  unit  would 
correspond  to  a  ring  unit. 

We  have  also  carried  out  calculations  for  the  monomer,  dimer,  trimer 
and  closed  6-ring  in  known  crystal  structures  (using  actual  coordinates 
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from  the  literature).  Table  V  lists  CE/q|  for  K  =  1  and  CE  for 
Q2  =  0.8  and  K  =  3.5  for  these  cases.  We  note  from  Eq.  (3)  that  CE 
varies  slowly  with  distance  and  expect  very  little  variation  within  the 
groups.  The  energy  differences  between  groups  are  noticeable;  i.e. 

-  12.5  Kcal/mole  between  monomer  and  dimer  or  between  dimer  and  closed 
6-ring.  The  energy  of  the  6- ring  can  be  considered  a  very  good 
approximation  for  the  energy  of  the  structure  as  a  whole  with  the  CE  of 
the  monomer  contributing  the  major  portion  of  this  energy.  Since  the 
standard  heat  of  formation  of  Si O2  from  the  elements  Si  and  0  ~  205  KPM, 
we  can  expect  the  covalent  energy  contribution  to  be  -  65  KPM  with  an 
Ey  -  70%  ionic.  If  we  had  reduced  Q2  to  0.7,  then  Ey  -  55%  ionic. 

With  either  of  these  choices,  AEy  for  the  closed  6-ring  between  the 
different  polymorphs  is  ~  2  -  4  KPM  which  is  the  experimentally  observed 
range.  Thus,  we  may  deduce  that  the  insight  provided  by  these  simple 
calculations  is  fully  valid  for  our  structural  assessment  and  we  may 
begin  to  consider  changes  associated  with  changing  the  important  angles 
0,  5p  and  in  the  dimer. 

The  dimer  can  be  defined  in  the  usual  way  via  0,  6^  and  5^  or  v''a 
the  alternate  scheme  described  in  an  earlier  report(^).  In  this 
alternate  but  equivalent  scheme,  three  oxygen  atoms  form  the  reference 
plane  (the  bridge  0  plus  one  0  from  each  tetrahedron).  This  scheme  is 
convenient  for  describing  ring  structures  in  terms  of  bridging  oxygens. 
The  dihedral  angles  (aj,  a2)  with  respect  to  this  reference  plane  and 
the  angle  O-OB-O  (<t>)  define  the  dimer  in  this  scheme.  We  have 
calculated  the  CE  for  various  values  of  three  angular  parameters  in  both 
schemes.  Some  of  these  results  are  summarized  in  Table  VI.  Four 
diagrams  have  been  presented,  one  for  each  of  <j>  =  210,  240,  270  and 
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300°.  In  the  horizontal  direction,  one  of  the  dihedral  angles,  ,  is 
varied  in  steps  of  30°  from  90°  to  270°.  In  the  vertical  direction,  the 
other  dihedral  angle,  ,  is  varied  from  0°  to  120°.  Any  intersection 
completely  defines  a  dimer  ($,aj  aj)  and  we  present  the  following 
calculated  information:  (1)  the  relative  CE  for  that  particular  dimer 
configuration,  (2)  the  minimum  0-0  separation,  rm,  between  the  two 
tetrahedra  in  the  dimer  and  (3)  the  corresponding  Si-O-Si  angle  (e)  . 

We  note  in  Table  VI  that  the  variations  in  E  =  K(CE)/7.81  Qjj  are 
generally  within  -  1  KPM  and  close  to  the  values  confuted  for  the 
dimers  in  known  structure  presented  in  Table  V.  The  shaded  areas  of 
Table  VI  are  the  important  ones  for  us  here.  They  correspond  to 
relatively  larger  values  of  CE,  small  values  of  rm  and  values  of 
0  3  130°  .  All  these  three  quantities  are  strongly  correlated!  This  is 

because,  as  the  tetrahedra  are  brought  closer  together,  the  0-0  distances 
are  decreased  and  9  decreases.  However,  because  the  CE  decays  only 
very  slowly  with  distance,  aCE  is  quite  small.  In  nature,  the 
separation  between  oxygens  from  neighboring  tetrahedra  are  seldom  less 
than  3.0  A  .  Even  under  very  high  pressure,  the  smallest  0-0  distance 
between  adjacent  tetrahedra  is  about  2.84  A  in  low  quartz.  Thus,  these 
shaded  areas  in  Table  VI  represent  dimer  states  with  a  very  low 
occurrence  probability. 

The  change  in  CE  due  exclusively  to  interactions  between  atoms  in 
different  tetrahedra,  aCE  ,  is  presented  in  Table  VII  for  some  randomly 
selected  values  from  the  dozens  of  dimer  configurations  analyzed 
according  to  the  (e.Aj.Aj)  scheme.  Because  of  the  weighting  scheme, 
the  monomer  contribution  in  a  dimer  should  be  taken  as  -96.21  KPM  (for 
K  =  1)  and  not  -81.35  as  given  in  Table  V.  Adding  the  aCE  =  5.36  (for 
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K  =  1)  from  Table  VII  due  to  interaction  in  the  dimer  gives  the  proper 
dimer  CE  (for  K  =  1)  of  -90.85,  which  is  in  close  agreement  with  the 
values  given  for  dimers  in  known  structures  (Table  V).  Once  again  the 
correlation  is  unmistakable  between  the  high  positive  values  of  aCE  and 
the  low  0-0  distances.  Very  low  values  of  0-0  would  mean  the  inter¬ 
section  of  tetrahedra  (tetrahedra  deformation  required).  It  is  also 
noted  that,  in  general,  a  aCE  >  6.00  for  (K  =  1)  corresponds  to 
e  Z  120°  while  a  aCE  <  5.00  for  (K  =  1)  corresponds  to  0  S  160°  . 

Although  it  seems  as  if  large  8  =  (0  ~  180°)  are  energetically 
favorable  for  the  dimer,  when  the  second  neighbor  interactions  are 
considered  (as  in  a  trimer)  the  energetically  favored  angles  revert  to 
the  range  of  9  ~  144°  .  The  non-bonding  interactions  will  also 

contribute  to  this  feature-favoring  a  closer  packed  structure  as  long  as 
the  separations  between  atoms  in  the  nearby  tetrahedra  are  not  less  than 
the  van  der  Waals  separation. 

Three  Membered  and  Four  Membered  Rings 

The  3-ring  of  oxygen  can  be  taken  to  be  an  equilateral  triangle 
($  =  60°)  .  Intuitively,  we  expect  the  Si  to  be  outside  the  triangle  so 
that  the  tetrahedra  do  not  intersect.  Taking  the  Si  of  one  tetrahedron 
to  lie  in  the  plane  of  the  0-ring,  a  few  combinations  of  the  dihedral 
angles  for  the  other  two  tetrahedra  were  investigated  and,  for  rm  5  3.0  , 
0  ~  125°  -  131°  generally  but  with  some  spread  down  to  0  -  110°  .  For 

the  4-ring,  we  assume  it  to  be  square  and  with  the  Si  of  one  tetrahedron 
lying  in  the  plane  of  the  ring.  The  remaining  three  tetrahedra  are 
defined  by  the  angular  parameters.  For  the  four  sets  of  chosen 
parameters  given  in  Table  VIII,  we  note  only  a  small  variation  in  CE  and 
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in  rm  while  120°  <  8  <  151°  .  From  this  we  can  deduce  that  a  given 
n-ring  will  exhibit  a  spectral  distribution  of  0  peaking  at  a  maximum 
probability  of  occurrence  for  a  minimum  Ey. 

A  more  detailed  account  of  the  study  of  trimers  will  illustrate  our 

approach  to  modeling  of  silica  much  better.  A  trimer  will  be  more 

indicative  of  the  ring  structure  that  is  likely  to  develop.  Three 
silicons,  Siy,  Si 2 .  and  Si 3  and  ten  oxygens  1-10  form  the  cluster. 

Oxygens  1  and  2  and  the  three  silicons  are  assigned  a  weight  =  1  and  the 
remaining  oxygen  a  weight  =  .5.  There  are  78  distances  in  all  among  the 
13  atoms  of  the  cluster.  Taking  the  central  unit  as  reference,  three 
angular  parameters  a,6,y  (Fig.  5)  are  chosen  to  define  each  of  the 

other  two  units.  In  other  words  we  define  two  dimers.  These  angles  are 

conveniently  chosen  as  the  Euler  angles.  The  choice  of  the  angular 
parameters  in  each  case  is  decided  by  the  ease  of  computation  and  the 
different  sets  of  parameters  can  be  interpreted  in  terms  of  one  another. 
The  angles  S  are  the  angles  between  the  sides  of  a  potential  ring 
structure. 

Systematic  tabulation  of  the  energy  of  the  trimer  as  a  function  of 
the  six  angles  defining  it  is  too  time  consuming,  may  be  difficult  to 
interpret  and  is  unnecessary.  Instead  we  decided  to  sample  the  structure 
over  the  entire  'space'  of  the  six  trimer  angles.  Random  values  in  the 
full  range  0-360°  were  assigned  for  the  six  trimer  angles.  Sieves  were 
used  so  that  the  full  calculation  (i.e.,  determining  the  78  distances  and 
finding  the  CE  for  the  trimer)  was  done  only  for  those  passing  through 
the  sieve.  For  these  initial  calculations  the  sieve  was:  when  any  0-0 
separation  was  found  to  be  less  than  a  chosen  minimum,  the  set  of  angles 
were  abandoned  and  a  new  set  of  random  values  chosen  and  the  procedure 
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repeated.  In  any  case  the  TOT  angles  were  calculated  from  Sij  -  Si 2  and 
Sij  -  Si 3.  For  the  initial  studies  we  chose  a  minimum  separation  between 
oxygens  =  2.8  A.  We  chose  this  value  because  no  known  polymorph  of 
silica  had  0-0  separation  less  than  this  value  and  because  of  the 
repulsive  energy  data,  using  the  Lennard-Jones  potential  parameters  of 
Table  IX,  indicate  a  major  energy  penalty  for  0-0  distances  less  than 
this.  This  choice  of  0-0  minimum  will  exclude  three  and  four  membered 
rings.  It  is  easy  to  recognize  that  at  least  one  of  the  9  distances  0 
(5-7)  to  0  (8-10)  should  be  zero  for  a  3-ring  and  about  2.6  A  for  a 
4-ring.  We  could  use  a  different  sieve  to  isolate  these  rings  in  the 
procedures  we  have  described.  In  the  test  runs  with  0-0  min  =  2.8  A,  we 
initially  tried  1200  sets  of  trimer  angles.  Of  them  only  61  passed  the 
sieve.  For  these  61  cases  all  the  78  distances  of  the  cluster  were 
calculated  and  the  CE  determined.  The  calculations  were  abandoned  in  the 
rest  of  the  tries  whenever  the  0-0  separation  fell  below  2.8  A.  The 
distribution  of  the  trimer  angles  over  the  1200  tries  was  examined  and 
found  to  be  uniformly  distributed  over  0-360°.  The  distribution  of  the 
TOT  angles  is  shown  in  Fig.  6.  The  range  of  the  TOT  angles  is  0-180°. 

The  TOT  angles  are  found  to  be  practically  uniformly  distributed  over  a 
large  range  20-140°.  It  must  be  recognized  that  the  TOT  angles  are 
functions  of  the  dimer  angles  and  uniform  distribution  of  the  dimer 
angles  over  the  range  0-360°  need  not  necessarily  imply  uniform 
distribution  of  the  TOT  angles.  However,  the  effect  of  the  filter  on  the 
distribution  of  TOT  can  be  seen  clearly.  The  TOT  angles  for  the 
successful  tries  are  distributed  over  a  much  smaller  range  110-180  with  a 
peak  around  140-150°. 


From  the  calculations  on  the  successful  tries,  regions  corresponding 
to  trimers  in  known  structures  can  be  identified.  For  instance,  in  a 
quick  scan  the  low  quartz  trimer  was  unequivocally  identified. 

The  continuation  of  this  work  will  take  the  following  course:  We 
can  use  N.B.E.  terms  (e.g.,  Lennard-Jones  6-12  potential)  instead  of  the 
coarse  sieve  of  minimum  0-0  separation.  Five,  6,7...  membered  rings  can 
be  built  using  similar  procedures.  This  is  not  a  diffiult  task.  For 
each  additional  unit  we  need  three  more  angles  to  be  specified.  We  can 
use  Euler  angles  as  we  did  for  the  trimer.  The  distribution  of  the 
3, 4, 5, 6, 7. . .membered  rings  can  be  studied. 

Alternately  we  can  increase  the  pool  of  trimer  angles  and  select  the 
angles  from  them  to  build  the  rings.  This  will  result  in  fewer 
rejections.  In  other  words,  this  pool  will  represent  an  initial  coarse 
filter.  From  any  of  the  angles  in  this  pool,  we  can  make  small 
variations  in  the  trimer  angles  and  find  the  energy  minima,  using 
optimization  techniques.  Programs  have  been  developed  by  us  for  these 
studies  while  working  on  a  different  molecular  simulation  problem. 
Modification  of  the  ring  structure,  e.g.,  introduction  of  new  species 
(vacancies  or  interstitials)  is  possible.  This  initial  work,  in  short, 
demonstrates  that  we  can  do  modeling  studies  on  silica  and  its 
modifications  using  relatively  simple  calculations  and  providing  useful 
information  on  the  energetics  of  formation. 

3.  Energetic  Aspects  of  Varied  Atomic  and  Molecular  Species 
in  Silica 

One  important  goal  of  using  a  molecular  orbital  approach  to  model 
the  energetic  aspects  of  interstitial  species  in  silica  is  determining 
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the  bridginq  oxygen  vacancy,  Vq,  and  Frenkel  defect,  Vq  +  Oj ,  formation 
energies.  This  information  is  essential  for  modeling  the  network  versus 
interstitial  diffusion  for  various  species  as  well  as  for  understanding 
the  mechanisms  of  formation  and  annealing  of  fixed  oxide  charge. 

The  use  and  limitations  of  the  CNDO/2  method  were  discussed  in  the 
previous  report^1).  The  most  serious  limitation  of  this  MO  technique  for 
modeling  vacancy  formation  is  that  bond  lengths  tend  to  be  overestimated 
for  second-row  elements.  CNOO/2  calculations  using  normal  bonding 
parameters  yield  equilibrium  Si-0  distances,  d(Si-O)  ,  greater  than  the 
values  found  in  natural  silicates.  For  example,  the  calculated 
equilibrium  Si-0  distance  using  normal  bonding  parameters  and  the  HgSi 2^7 
cluster  is  greater  than  2.0  A  (3),  Pereas  measured  values  for  silica 
range  from  1.6  to  1.63  A  .  Thus,  when  the  bridging  oxygen  is  moved 
from  its  normal  position  in  the  HgSi^Og  cluster,  the  total  energy 
decreases  as  d)Si-O(br))  increases.  This  leads  to  unreasonably  low 
total  energies  for  a  HgSi 2^7  cluster  without  a  bridging  oxygen,  Vq  . 

To  decrease  the  calculated  equilibrium  values  of  d(Si-O)  ,  it  was 
necessary  to  increase  the  orbital  exponent  of  Si,  £<j-j  >  ^rom  the  normal 
value  of  1.383(3)  to  1.86.  Increasing  the  bonding  parameter, 

^ ,  contracts  the  Si  orbitals  and  produces  a  minimum  in  the  total 
cluster  energy  at  d(Si-O)  =  1.61  for  =1.86  .  Although  the  goal  of 
matching  calculated  and  measured  equilibrium  d(Si-O)  values  using 
CNOO/2  was  attained  by  modifying  ,  the  electronic  structure  of  the 

molecule  was  altered  as  well.  For  example,  the  charge  on  the  bridging 
oxygen  became  less  negative  (from  -0.75  to  -0.52)  indicating  an  increase 
in  the  covalent  character  of  the  Si-0  bond.  High  bridging  oxygen  atomic 
charges  were  also  found  by  Newton  and  Gibbs^)  using  ab  initio  methods  to 
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study  HgSi 2O7  clusters  and  including  d-orbitals  in  the  basis  set.  No 
molecular  orbital  method  applied  to  HgSi 2O7  appears  to  simultaneously 
yield  values  of  atomic  charges  and  bond  lengths  completely  consistent 
with  the  experimental  values  for  silica. 

With  the  Si-0  bond  length-Ey  relationship  optimized,  it  was  possible 
to  calculate  energies  of  formation  of  vacancy  species.  Bridging  oxygen 
vacancy  formation  in  a  Si-0  network,  Vq,  is  modeled  by  the  following 
reaction: 


H6Si20;  t  HgSjOgVg  +  0  ;  AEy  =  +  8.967  eV  .  (5) 

The  HgSi207  cluster  models  a  Q(br)  site  in  an  Si-0  network 
whereas  HgSi20gV0  represents  the  vacancy  configuration,  V0.  For 
reaction  (5),  the  equilibrium  Si-0  bond  distance  was  calculated  to  be 
1.61  A  at  =  1.86  .  An  eclipsed  dimer  configuration  similar  to  that 

discussed  in  the  previous  report  was  used  fixing  <  Si-O-Si  *  143.6° 
and  <  Si-O-H  =  180°  .  The  vacancy  configuration  consists  of  the  same 

dimer  configuration  without  a  bridging  oxygen.  The  high  positive 
reaction  energy  for  (5)  suggests  a  low  vacancy  density  would  be  produced 
by  a  reaction  involving  only  neutral  species  and  sites  and  leading  to  the 
release  of  monatomic  oxygen.  Performing  a  similar  calculation  for  the 
tetramer  species  using  ^ .  =  1.86,  d(Si-O)  =  1.62  and  d(O-H)  *  0.96 
led  to  aEt  =  ♦  8.8984  eV  .  This  is  a  remarkable  correspondence  with 
the  results  of  Eq.  (5)  considering  the  differences  in  d(Si-O)  and  d(0- 
H)  between  the  dimer  and  tetramer. 

A  much  more  favorable  reaction  for  the  formation  of  0(br)  vacancies 
is  that  involving  the  reaction  product  02;  i.e.. 
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2Hg Si 2 0?  i  2HgSi 20g VQ  +  02  ;  aEt  =  +0.4972  eV  (6) 

where  the  equilibrium  d(0-0)  value  in  O2  gas  is  1.13  A  at  0°K.  In 
this  case  the  energy  of  formation  at  0°K  for  Vg  at  the  free  Si 02  surface 
is  only  0.25  eV.  If  we  assume  that  aEj(T)  =  aEj(0°K)  and  use  the 
equilibrium  relationship 

a£t(T)  =  -RT  *n[(Xv  /X0(bp))2Pg  ]  =  11.466  kcal/mole  (7) 

We  find  that  about  6%  of  the  bridging  oxygen  sites  at  the  free  surface 
would  be  vacant  at  1000°K  for  Pg  =  1  atm  .  Since  the  energy  to 
form  a  vacancy  in  a  tetramer  was  found  to  be  almost  identical  with  that 
for  a  dimer  and,  since  in  our  previous  report  we  found  that 
aEj  =  +  7.387  eV  (Eq.  (8))  for  placing  an  02  in  the  center  of  a 
tetramer,  we  expect  the  following  to  hold 

3  tetramer  £  2(tetramer-Vg)  +  (tetramer  with  02  in  center);  aEj  »  +7.884  eV 

(8) 

so  we  have  for  the  mole  fraction  Xv  of  bridging  oxygen  vacancies  in 

v0 

the  bulk  of  the  tetramer  silica 

Xv  .  exp(-60/RT)  (9) 

v0 

with  RT  in  kcals.  This  would  lead  to  a  fraction  of  about  10"^  vacant 
bridging  oxygen  sites  in  the  bulk.  Since  we  expect  the  adsorption 
energy  of  an  02  molecule  at  the  exterior  (obtuse  angle)  of  a  tetramer  to 
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be  much  less  than  in  the  interior  (acute  angle)  of  a  tetramer,  the 
actual  formation  energy  for  a  bridging  oxygen  vacancy  in  bulk  silica  (6- 
ring  rather  than  4-ring)  will  be  appreciably  smaller  than  given  by  Eq. 
(5). 

In  Table  X,  using  £  =  1.86  ,  the  energy  change  associated  with 
moving  an  oxygen  from  its  bridging  oxygen  site  in  a  tetramer  to  the 
center  of  the  ring  is  given.  Thus,  for  a  Vg  +  Qj  in  the  same  ring, 
the  energy  of  formation  is  5.94  eV.  However,  by  using  Eq.  (4)  plus  the 
energy  change  due  to  O2  dissociation  in  a  tetramer;  i.e.. 


tetramer  +  tetramer  with  0^  at  (0,0)  *  2(tetramer  +  0  at  (0,0));  AE-p  =  4.56 

(10) 

from  our  previous  report  (1)  (Eqs.  (8)  and  (9)),  the  energy  for  forming  a 
dissociated  Frenkel  defect  in  a  tetramer  solid  is  given  as 

4  tetramers  +  2(tetramers  -  Vg)  +  2(tetramers  -0j);  AE-p  =  +3.32  eV 

(ID 

so  that  the  energy  to  form  a  dissociated  Frenkel  defect  is  +1.61  eV  and 
this  type  of  dissociated  defect  is  stable  in  such  a  tetramer  solid. 

In  the  calculations  involved  with  removing  a  bridging  oxygen,  its 
charge  changed  from  ~  -0.5  to  -0.2  at  1  A  displacement  and  to 
-  -0.015  at  1.9  A  displacement  so  that  the  charge  on  the  removed  0  and 
upon  the  remaining  fragment  shrink  to  zero  as  the  0  becomes  completely 
separated.  The  reaction  of  +2  charged  bridging  oxygen  forming  a  charged 
vacancy  is  much  more  favorable  energetically  than  Eq.  (5);  i.e., 
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(12) 


*  Iv 


H6S12°7+  *  H6Si2°6V0+  +  0  :  aEt  =  +  0,279  eV 

However,  considering  the  formation  of  a  +2  charged  0(br)  site  including 
the  electron  work  function 

H6Si2°7  *  H6Si2°7+  +  2e"  ’  aET  =  +  45-586  eV  (13) 

This  tells  us  that  Eq.  (12)  is  not  favored  because  the  initial  charged 
species  is  highly  unfavored. 

Frenkel  defect  formation  energies  have  also  been  calculated  for  both 
neutral  and  charged  configurations;  i.e., 

2HgSi 20?  X  HgSi^gOj  +  rtgSi 2°6 Vo  ’  aET  =  '  1,065  eV  * 

(14) 

Reaction  (14)  models  the  formation  of  an  oxygen  interstitial  ( Hg Si ) 0 j 
and  oxygen  vacancy  (HgSigOg)!/^  simultaneously  from  two  bridging  oxygen 
sites.  The  Ej  of  the  oxygen  interstitial  configuration  was  calculated 
using  the  HgSi^Oy  configuration  as  above  with  an  oxygen  adatom  at  an 
equilibrium  distance  of  1.25  A  from  the  bridging  oxygen  along  z  , 


OH 

OH 


0  0 

H  H 


HgSljOg  =  HgSigOyOj 
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The  negative  reaction  energy  of  Eq.  (14)  suggests  that  the  large 
negative  energy  of  formation  of  0  interstitials  more  than  compensates 
for  the  unfavorably  positive  reaction  energy  for  vacancy  formation  and 
that  such  uncharged  Frenkel  defects  will  form  spontaneously  at  the  free 
surface.  In  the  interior  of  the  volume,  for  6-rings,  the  formation  of 
these  neutrally  charged  species  is  expected  to  involve  a  AEy  >  0  . 

For  charged  Frenkel  defects,  we  have  two  important  possibilities 


2H6Si20?  ♦  H6Si  2°7°I~  +  H6Si2°6V0+  ;  =  +  21  *75  eV 

(15) 

and 

4HgSi207  J  2H6Si  2070j"  +  HgSi^OgVg  +  HgSi20gV^+  +  2e"  ;  aEt  =  +  16.09  eV 

(16) 

In  the  case  of  reaction  (16)  where  two  Frenkel  defects  are  formed,  the 

oxygen  vacancies  of  zero  charge  and  +2  charge  might  be  used  to  model  the 

formation  of  two  vacancies  each  with  a  charge  of  +1  leading  to  a  Frenkel 

defect  formation  energy  ~  +  8  eV  for  this  defect.  Clearly,  these 

charged  defects  do  not  appear  to  be  energetically  favorable  and  their 

presence  will  be  enhanced  only  at  the  higher  temperatures,  low  Pn  and 

UZ 

in  the  presence  of  reactions  that  consume  electrons. 

Since  the  Ey  of  HgSi 2O7  could  be  minimized  at  d(Si-O)  =  1.61  by 
changing  the  bonding  parameter,  £  ,  the  interaction  energies  determined 
previously  were  recalculated  to  ascertain  if  any  differences  in  reaction 
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behavior  would  result  from  the  new  parameterization.  It  was  found  that 

the  Ey  j/s_.  distance  relationships  for  02  interactions  with  0(br)  in 

HgSi 2O7  are  similar  for  both  normal  and  altered  parameterizations.  For 

example,  in  the  configuration  with  the  0-0  axis  perpendicular  to  the  Si-0- 

Si  plane,  the  equilibrium  0(br)-admolecule  distance  is  1.15  A  for 

C  =  1.86  and  1.2  for  the  normal  parameter.  However,  for  the  interaction 

of  Si  atoms  with  bridging  oxygens,  much  larger  differences  were  observed 

in  Ey  _vs_.  equilibrium  distance  relations  between  the  two 

parameterizations.  For  neutral  HgSijOy  ,  the  equilibrium  Si(adatom)- 

0(br)  distance  was  1.76  A  for  5  =  1.86  but  greater  than  2.0  A  for  the 

2+ 

normal  parameterization.  Likewise  d( Si -0(br) )  for  H^Si^Oy  was 

1.58  A  using  the  new  orbital  exponent  but  also  greater  than  2  A  for  the 

normal  value  of  5  . 

The  results  of  0(br)  adatom/admolecule  interactions  suggest  that 
changing  will  only  affect  the  interaction  behavior  between  0(br)  and 

Si  and  not  between  0(br)  and  other  admolecules.  This  result  also  confirms 
the  importance  of  bonding  interactions  between  bridging  oxygens  and 
adatoms/admolecules  as  opposed  to  interactions  between  adatoms/admolecules 
and  other  atoms  such  as  Si  in  the  HgSi20y  cluster.  This  conclusion  is 
significant  since  it  supports  the  validity  of  CND0/2  calculations  using 
normal  Si  bonding  parameters  for  modeling  interactions  between  bridging 
oxygens  and  adatoms  in  aqueous,  gas,  and  crystalline  environments. 

Further  adatom/admolecule  -  0(br)  interaction  studies  are  underway  to 
verify  and  extend  the  above  results  and  conclusions.  The  admolecules 
include  Ct.,,  N2>  02  H2  F2>  H20  and  C02  using  both  parameterizations. 

The  results  for  the  normal  parameterization  (c  =  1.383)  are: 
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h  H6S12°7  5 

H6Si207CA2(ll); 

AEy  =  -9.509  kcal/mole 

(17) 

n2  ♦ 

H6Sl  2°7  + 

H6Si207N2(ll)  ; 

aEt  =  -6.627  kcal/mole 

(18) 

°2  + 

H6Sl  2°7  > 

HgSi20702(ll)  ; 

aE^  =  -2.083  kcal/mole 

(19) 

H2  + 

H6Si2°7  + 

HgSi207H2(ll)  ; 

aEj  =  -0.603  kcal/mole 

(20) 

F2  + 

H6Si2°7  + 

H6Si2°7F2(11)  ; 

aEj  =  -0.282  kcal/mole 

(21) 

The  equilibrium  distances  between  the  diatomic  gas  axis  and  0(br)  for  the 
molecules  on  the  right  of  Eqs.  (17)- (24)  are,  respectively  2.4  A  , 

2.1  A.  2.1  A,  2.5  A  and  2.7  A  .  For  contrast,  the  equilibrium  distances 
in  the  diatomic  gases  at  8°K  are  d(N-N)  =  1.14  A,  d(Cz-Ci)  =  1.98A, 
d(0-0)  =  1.13  A,  d(H-H)  =  0.746  A  and  d(F-F)  =  1.12  A.  .  The  minimum 
energy  parallel  configuration  has  the  axis  of  the  diatomic  molecule 
parallel  to  the  Si-O-Si  plane. 

Table  XI  illustrates  the  relationship  between  the  silicon-oxygen 
bridge  bond  energy,  E$i-0(br)’  and  AEy  for  the  reactions  17-21  as  well 
as  other  admolecule  reactions.  These  results  indicate  that  Vq  formation 
may  be  favored  by  reactions  with  the  most  negative  AEy  since  the  Si- 
0(br)  bond  energy  increases  significanty  as  aEj  decreases. 
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TABLE  I 


Argon 

Ion-Etching 

Data  for 

Oxide  Removal 

from  Sil icon 

Distance 

(inches) 

Vol tage 
(kV) 

Argon 

(pm) 

Etch 

Time 

(min) 

Thickness 

Change 

o 

(A) 

Etch 

Rate 

0 

(A/min) 

4 

0.75 

4 

30 

no  ie 

- 

4 

0.75 

8 

40 

none 

- 

4 

2.0 

4 

20 

245 

12.25 

4 

2.0 

4 

40 

650 

16.2 

4 

2.0 

8 

20 

270 

13.5 

4 

2.0 

8 

40 

580 

14.5 

5 

0.75 

4 

20 

none 

- 

5 

0.75 

4 

30 

none 

- 

5 

0.75 

4 

40 

none 

- 

5 

0.75 

8 

20 

none 

- 

5 

0.75 

8 

40 

none 

- 

5 

2.0 

7.5 

15 

none 

- 

5 

2.0 

7.5 

20 

none 

- 

5 

2.0 

7.25 

40 

130 

3.25 
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Data  on  Reactively  Sputtered  Si09  Films  {Ion  etched  substrate)  -  Substrate  Temperature  Dependence 
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TABLE  IV 


Relative  X-ray  Peak  Height  Data  for  Various  SiC  Films 


Film# 

Film/Substrate 

Process^ 

Relative  X-ray 

Peak  Height 

177 

SiC/(100)Si 

600W,  720°C,  3" 

1 .0 

178 

SiC/(lll)Si 

600W,  660° C,  4" 

5.3 

180 

SiC/tyco 

6Q0W ,  660°C,  4" 

0.84 

184 

Si C/ (111 )Si 

200W,  700°C,  4" 

1.0 

210 

SiC/(100)Si 

600W,  25°C,  4" 

0.0 

★ 

Tyco  ribbon  sapphire 

4. 

SiC  Source  power,  substrate  temperature,  source-substrate  distance 


TABLE  V 


Coulomb  Energy  for  Structural  Elements 
in  Different  Si02  Polymorphs 


El ement 

ce/q22  (KPM) 

(K=  1) 

CE(KPM) 

(K  =  3.5) 

02=0.8 

Monomer  . 

Regular  Tet 

-639.5 

-117.05 

Low  Quartz  (105PA) 

-626.6 

-115.72 

"  "  (68  PA) 

-633.8 

-116.00 

"  Crystobal ite  (28°C) 

-636.5 

-116.50 

"  "  (230°C) 

-639.5 

-117.05 

High  "  (300°C) 

-634.5 

-116.14 

”  "  ( 1 200°C) 

-636.3 

-116.47 

Dimer 

Low  Quartz  (10b  PA) 

-701.3 

-128.36 

"  "  (68  PA) 

-701.2 

-128.34 

"  Crystobal ite  (28°C) 

-707.0 

-129.40 

High  "  (300°C) 

-704.8 

-129.00 

Trimer 

Low  Quartz  (68  PA) 

-723.7 

-132.46 

Closed  6-Rinq 

Low  Quartz  (105  PA) 

-769.7 

-140.88 

"  (68  PA) 

-769.2 

-140.80 

Low  Crystobalite  (28°C) 

-776.6 

-142.16 

"  "  (230°C) 

1 

CO 

—4 

o 

-142.96 

High  "  (30Q°C) 

-774.5 

-141 .77 

"  "  ( 1 200°C) 

-776.0 

-142.04 

High 


TABLE  VII 


ACE  Contributions  Due  to  Interactions 
Between  Adjacent  Tetrahedra 


Random  ACE  (KPM) 

Sample  (for  K=l)  ^  (KPM)  (0o‘0)min 


TABLE  VIII 

Variation  of  9  .  ,  r  and  CE  for  the  4-Ring 

mi  n  m 


CE  (for  K=1)  KCE 

Angles  Q2=0.8  ~q|~  (0;  0)min  9min 

(Deg)  (KPM)  (KPM)  A  (Deg) 

120  60  180  -99.0  -774  3.47  135 

150  30  180  -99.9  -780  3.70  151 

120  300  150  -99.0  -774  3.47  120 

150  330  180  -99.9  -780  3.70  151 


TABLE  X 


Undissociated  Frenkel  Defect  Calculation  (£  -  1.86) 


2  *  1.87 


\  / 


1 

Ey(a.u. ) 

187 

-245.208426 

1.7 

-245.198887 

1.5 

-245.192528 

1 .3 

-245.189946 

l  .0 

-245.165820 

0.7 

-245.090046 

0.3 

-245.013903 

0.0 

-244.984417 
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TABLE  XI 


Admolecule-  Dimer  Interaction  Energies  for  Silica 
Adspecies  AET(Kcal/mole)  ESi-0(br)^eV^ 


F2 

-0.3 

-17.071 

H2 

-0.7 

-17.019 

°2 

-2.1 

-16.565 

N2 

-6.6 

-  16.286 

h2o 

-9.0+ 

-16.796 

-9.5 

-16.886 

OH' 

-17. 0+ 

-12.946 

HQ 

-28.7 

-14.565 

* 

4- 

-14.256 

C02(I) 

-35 

★ 

+ 

C02( II) 

-4.6 

-12.480 

H+ 

-317+ 

-9.308 

*,+  Configuration  and 

( 5 ) 

and  Brown v  ' 

energy  respectively  described 

in  deJong 

«  __  ^  i 

k  '  '  h  rijhaaiiiiiiSifiBBBiiaJi 
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FIGURE  CAPTIONS 


(a)  Arrangement  of  deposition  source,  Ar-cleaner  ,  oxygen  supply  ring  and 
substrate  heater  used  in  the  reactive  sputtering  of  Si 0X  films. 

(b)  Expanded  view  of  the  quartz  ring  for  supplying  02  to  the  rectangular 
substrate  surface. 

Si 0  film  thickness  at  constant  Pn  as  a  function  of  substrate  tempera- 

X 

ture  for  fixed  Si  source  power. 

Single-source  SiC  sputtering  system  with  Ar  ion-gun  for  substrate  sur¬ 
face  cleaning,  and  with  "portable"  tin  source  situated  in  most  likely 
position. 

Electron  microprobe  results  of  one  dual-source  and  three  single-source 
SiC  depositions  showing  atomic  percentages  of  Si  and  C  as  a  function  of 
sample  position. 

(a)  Schematic  representation  of  the  three  tetrahedra  involved  in  a  silica 
trimer  (3  Si  +  IOC). 

(b)  Illustration  of  possible  angular  variations  a,  B,  y  for  the  three 
tetrahedra  with  one  held  fixed  (6  degrees  of  freedom). 

Plot  of  the  number  of  trimers  as  a  function  of  bridge-bond  angle  0.  The 
upper  histogram  is  for  all  the  choices  in  the  sample  (total  =  1200); 
the  lower  histogram  is  for  the  choices  passing  through  the  energy  sieve 
dQ_0  7  2.8  A0  (total  =  61 ). 
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FIGURE  2 


-  43  - 


ion 

gun 


FIGURE  3(a) 


44 


rototing  toble 


atomic  percent  sili 


T 


0  12  3  4 

cm  along  sample 


FIGURE  4 


